Rabi oscillations, coherent properties and model qubits 
in two-level donor systems under terahertz radiation 
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Quantum confinement, magnetic-field effects, and laser coupling with the two low-lying states of 
electrons bound to donor impurities in semiconductors may be used to coherently manipulate the 
two-level donor system in order to establish the appropriate operational conditions of basic quan- 
tum bits (qubits) for solid-state based quantum computers. Here we present a detailed theoretical 
calculation of the damped Rabi oscillations and time evolution of the Is and 2p+ donor states in 
bulk GaAs under an external magnetic field and in the presence of terahertz laser radiation, and 
their influence on the measured photocurrent. We also detail the possible experimental conditions 
under which decoherence is weak and qubit operations are efficiently controlled. 



PACS: 42.65.-k, 71.55.Eq and 73.20D.Hb 

Some two decades ago, it was suggested that quantum- 
mechanical computers (QC) might provide much more 
computing power than machines based on classical 
physics 1 . More recently, the discovery of fast quan- 
tum algorithms for tasks such as prime factorization and 
searching disordered databases 2 as well as quantum er- 
ror correction codes 3 have arisen considerable interest in 
the search for suitable physical systems on which reli- 
able quantum hardware could be realized. A precisely 
defined two-fold Hilbert space may be regarded as the 
base of a QC, under the conditions that decoherence is 
weak and single-qubit and two-qubit unitary operations 
are controlled. 

Many implementations have been suggested for quan- 
tum computation involving manipulation of nuclear 
spin states in bulk solution using nuclear magnetic 
resonance techniques (NMR) 4 , trapped ions 5 , photons 
in microcavity 6 , and specifically selected donor 7 or 
exciton 8-10 states in semiconductor systems as qubits to 
be manipulated by laser radiation. Solid-state propos- 
als open up the possibility of fabricating large integrated 
networks that would be required for realistic applications 
of quantum computers. Two interesting suggestions for 
such semiconductor QCs are based on qubits using the 
spin of a phosporus donor nucleus in bulk silicon 11 and 
the spins of electrons trapped in a GaAs quantum dot 
(QD) 12 ' 13 . Of course, a silicon-based computer is ex- 
tremely attractive from many aspects but unfortunately 
difficulties concerning its practical realization have been 
pointed out recently 14 . The main problem is related to 
the strong oscillations in the exchange splitting- basic in- 
gredient for the two-qubit operation in silicon-based QC, 
originated from degeneracy of the conduction band mini- 
mum in Si. This situation does not occur in GaAs where 
coupled-QD structures have been previously theoretically 
investigated 12,13,15 considering two coupled GaAs QDs. 
This corresponds essentially to "a 2D - hydrogen-like 
molecule" , in the presence of an external magnetic field. 



An extensive study by Hu and Das Sarma 13 and de Sousa 
et al 15 suggested that the two-QD system would provide 
the necessary two-qubit entanglement required for QCs. 
In this case, the QDs provide the tag for each qubit which 
are represented by the spin of the QD-trapped electron, 
and single-qubit operations are performed on the tag-spin 
by a local external magnetic field or equivalently by an 
efficient manipulation of the local spin states 16 . The two- 
qubit operation, which rotates one of the tag-spin qubit 
in a precise angle if, and only if, the other control-qubit is 
oriented in a well defined unchanged direction, is realized 
via exchange interaction between the electrons in each 
QD. This means that in order to have a precise control 
of two-qubit operations it is necessary to fine-tune the ex- 
change coupling through gates which are physically real- 
ized by external fields. One should note that a QC based 
on this architeture would be effective only if the spin de- 
coherence time is much longer than the time involved in 
the single- and two-qubit operations, which could pose 
a problem, as these operations would be controlled by 
switching electric and/or magnetic fields and this can- 
not be performed very fast. It is possible to overcome 
this limitation if, instead of slow external gate potential 
and/or magnetic fields, one uses a laser-probe pulse to 
control the qubit operations. This can be achieved via 
coherent manipulation of two-level systems through Rabi 
oscillations in donor 7 or exciton 8-10 states, by applying 
electromagnetic fields. Cole et aU and Zrenner et al 10 
have demonstrated that the coherent optical excitations 
in a QD and bulk GaAs two-level systems, respectively, 
can be converted into deterministic photocurrents. Here, 
the impurities (dots) provide the tag for each qubit which 
are represented by the two-level system. Moreover, it has 
been shown 17 that donor states in QDs, under magnetic 
and pump-laser 18 fields, may be properly manipulated in 
order to produce the required robust donor states for ter- 
ahertz coherent manipulation of qubits. In the present 
work, we present a detailed theoretical calculation of the 
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time evolution of the Is and 2p+ donor states in bulk 
GaAs and their influence on the photocurrent 7 , and thor- 
oughly detail the conditions under which decoherence is 
weak. 

A set of optical Bloch equations 18 is used to describe 
the time evolution of the elements of the density matrix 
within a two- level model for the impurity system, i.e., 



dt 

dp22 

dt 

dt 
dp2i 
dt 
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where is the THz laser frequency, W21 is the energy 
separation of the Is and 2p + impurity levels, Qr = 
ETHzdf 2 /h is the Rabi frequency, Ethz is the ampli- 
tude of the terahertz electric field (in the z-direction), 
and df 2 is the ^-component of the ls-2p + dipole matrix 
element 19 . The parameters 71, 72, and 73 are recombina- 
tion rates as defined in Cole et aP . The Is- and 2p + -like 
impurity states are calculated by using hydrogenic-likc 
wave functions with exponentially-decaying variational 
parameters taken as Ai s and A 2p+ . The calculated ener- 
gies, in effective units, are given by 



ei s = E c 



2 



72- - — + o7 A ls 



£ 2p+ 



\ 2 

A 2p+ 



A 2 



p+ 



+ 7+^7%+ 



(2) 



(3) 



where E c = E° + 7, E c being the magnetic- field depen- 
dent conduction-band edge, and 7 = ehB/2m*cR* is the 
effective magnetic energy, with R* w 5.9 meV the GaAs 
effective Rydberg. The parameters Ai s and \2p + are ob- 
tained variationally as a function of the applied magnetic 
field (in the z-direction). Fig. 1(a) shows the energies 
of both impurity states together with the magnetic-field 
dependence of the energy of the conduction-band edge. 
By using the value of v = 2.52 THz for the radiation 
frequency as in the experiment by Cole et aP , we cal- 
culated the magnetic field value of 3.4 T for which the 
ls-2p + transition is tuned with the radiation. This is in 
good agreement with the experimental value of 3.51 T 
obtained by Cole et aP at very low values of the tera- 
hertz electric field. The Is and 2p + hydrogenic-like wave 
functions, for the impurity under a magnetic field of 3.4 
T, are used to theoretically evaluate the Rabi frequency 
as a function of the THz electric field, and results are 
shown as a full curve in Fig. 1(b). Also shown are the fit 
parameters by Cole et aP as full triangles. We note that 
Cole et aP have pointed out that the uncertainty of the 
absolute scale of Ethz is ± 50 % and the uncertainty in 
relatives values of Ethz is negligible. With this in mind, 



one observes that the open triangles, which correspond 
to a horizontal shift of the fit parameters by 1.1 x 10 4 
Vm _1 , are in quite good agreement with our theoretical 
prediction for the Rabi frequency. 

The equations of motion of the density matrix may 
be solved via the rotating- wave scheme 18 , and solutions 
essentially oscillate with a frequency close to the Rabi 
frequency, with damping terms depending on the various 
recombination mechanisms. To discuss the photocurrent 
measurements in Fig. 4(a) of Cole et aP , we choose 71 =0, 
72 = 1.5 x 10 11 rad s" 1 , 73 = 1.0 x 10 11 rad s" 1 , and 
Qr = 3.0 x 10 11 rad s _1 . Our calculated results for 
the photocurrent [taken as proportional to 1 — pn(t)] in 
Fig. 2(a) show that this choice of parameters fits quite 
well their experimental data for the photocurrent, in the 
range of measured detunings. Note that the theoreti- 
cal photocurrent results in Fig. 2(a) indicate [cf. Fig. 
1(b)] that the experimental data in Fig. 4(a) of Cole et 
aP should be assigned to a laser amplitude of Ethz ~ 
3.1 x 10 4 Vm" 1 . In Fig. 2(b), we display the theo- 
retical photocurrent for a laser detuning S— -0.014 THz 
as compared with p22(t) (number of donors at the 2p + 
state), and pn(t) + P22{t) (total population in the Is 
and 2p + states) . One notes that the populations of both 
the Is and 2p + states have a fast decrease with pulse 
duration. Therefore, it is clear that the measured pho- 
tocurrent may have a dominant contribution from the 
charge ionized from the 2p + state into the conduction 
band, making the existence of coherent Rabi oscillations 
more difficult to be observed, i.e., the fact that the 2p + 
state is immersed in the conduction band results in a 
strong decoherence process. It is important to mention 
here that it is possible to increase the confinement of 
the donor electrons by trapping the impurities in QDs 
or via the use of a pump-laser tuned below the semicon- 
ductor gap, in order to lower the energy of the 2p + state 
and to bring it below the conduction-band continuum 17 . 
This situation corresponds to a choice of the 73 ioniza- 
tion rate equal to zero, as shown for the calculated pho- 
tocurrent in Fig. 3(a). Alternatively, one may efficiently 
control the decoherence process through a diluted dop- 
ing procedure in which the inhomogcncous broadening 
of the ls-2p + transition energy is conveniently reduced. 
This occurs because the overlap between the impurity 
wave functions diminishes. Again, the confinement of 
the electronic states via a pump-laser field may be used 
to mimic this effect. Such behaviour may be modelled 
by assuming the total dephasing rate 72 = 0, as dis- 
played in Fig. 3(b). One should note that a too low 
impurity concentration may result in a weak photosignal 
which is a rather undesirable experimental situation. In 
the theoretical photocurrent displayed in Fig. 3(c), it is 
quite clear the obvious fundamental role played by the 
strength of the electromagnetic field, through the Rabi 
oscillations, on the coherence process. 

In conclusion, we have obtained the solutions of the op- 
tical Bloch equations within the rotating-wave approxi- 
mation in order to discuss the proposal of Cole et aP . The 
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donor Is and 2p + hydrogenic-like states, in the presence 
of a homogeneous magnetic field, are obtained through 
a variational procedure. Theoretical calculations have 
shown that the resonance condition - for which the ls- 
2p + transition is tuned with the terahertz radiation - is 
achieved for a magnetic field in good agreement with the 
measured experimental value. Moreover, the theoretical 
results obtained for the THz electric field-dependence of 
the Rabi frequencies are found in excellent agreement 
with the experimentally fitted curve, provided that the 
uncertainty of the absolute experimental scale for the 
electric field is taken into account. The theoretical calcu- 
lations concerning the time-dependent behaviour of the 
photocurrcnt, a fundamental observable for an efficient 
determination of qubit operations, show that it is domi- 
nated by the ionized charge in the conduction band, and 
that this effect relates to the 2p + state being immersed in 
the conduction band. In that respect, we have detailed 
the possible experimental conditions for which decoher- 
ence effects are minimized so that qubit operations are 
efficiently controlled. 
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FIG. 1. (a) Theoretical magnetic - field dependence of the 
Is and 2p + donor states together with the field-dependence of 
the GaAs E c conduction-band edge. At 3.4 T, one obtains a 
ls-2p+ transition tunable to an energy corresponding to 2.52 
THz; (b) Rabi frequency as a function of the THz electric 
field (full curve), obtained via a variational calculation (see 
text). Full triangles are the fit parameters by Cole et al 7 , with 
open triangles corresponding to a horizontal shift by 1.1 x 10 4 
Vm" 1 . 



FIG. 2. (a) Theoretical 1 — pn(t) photocurrent (arbitrary 
units) as a function of pulse duration for various laser de- 
tunings <5=0 THz, -0.014 THz, -0.029 THz, and -0.043 THz; 
(b) Photocurrent for a laser detuning S— -0.014 THz as com- 
pared with pii(£) (number of donors at the 2p_|_ state), and 
pn(i) + P22{t). Results are calculated using 71 =0, 72 = 1.5 x 
10 11 rad s~\ 73 = 1.0 x 10 11 rad s _1 , and D. R = 3.0 x 10 11 
rad s _1 . 
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FIG. 3. Dependence of the photocurrent (arbitrary units) 
on the pulse duration for a fixed laser detuning <5 = -0.014 THz 
and (a) 73 = 0, and 73 = 1.0 (with 72 = 1.5 and £Ir — 3.0 ), 

(b) 72 = 0, and 72 = 1.5 (with 73 = 1.0 and Qr — 3.0); and 

(c) n n . = 3.0, and fL n = 4.0 (with 72 = 1.5 and 73 = 1.0). 
Rabi frequencies and rate parameters are given in units of 
10 11 rad s _1 . 
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